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Abstract 

A new mechanism for the primary photoinduced charge separation in photosynthesis 
is proposed. It involves as real intermediate between the excited special pair state P* 
and the primary charge separated state P^HJ^ a trip-trip-singlet P^Bj^, which consists 
of a triplet on the dimer P and a further triplet on the monomer Bl. Both combine 
to a singlet. The electron transfer is caused by spin exchange couplings. The transient 
spectrum of the short lived intermediate, formerly taken as evidence for the charge 
transfer state P~^B]^, is reinterpreted as a transient excitation of this trip-trip singlet. 



1 Introduction 

Photosynthetic systems convert photon energy into electrostatic energy within a few 
picoseconds with great efficiency. The initial charge separation takes place in a reaction 
center. Its basic structure (Deisenhofer et al 1984) (Fig. 1) seems to be preserved 
among all photosynthetic systems. Six pigments are arranged in an approximate C2 
symmetry forming an L and an M branch. Two make the special pair dimer. Two 
referred to as monomers are linked to it within van der Waals contact followed by 
two further pigments both being close to one of the monomers. The one on the L 
branch acts as the initial acceptor. In the special case of Rps. viridis, on which we 
want to base our calculations, the dimer P consists of two bacteriochlorophylls Pl 
and Pm noncovalently bound via 7r-orbital interactions of their pyrrol rings I. The 
monomers are two bacteriochlorophylls denoted Bl and Bm ■ They point towards 
the rings I of the dimer P with their rings III. The rings I of Bl and Bm face the 
corresponding rings I of the bacteriopheophytines Hl and Hm respectively. We will see 
that these basic structural features are essential to bring the energy location of the trip- 
trip singlet P^BJ^ below that of the initially excited dimer state P*. Furthermore we 
show that the spin exchange couplings are particularly favorable for this arrangement. 
This is not so for the commonly postulated intermediate charge transfer state P^Bl 
(Fischer and Scherer 1987; Scherer and Fischer 1989a; Scherer and Fischer 1989b; 



Holzapfcl ct al 1990; Scherer 1989). We predict it to lie above P* by about 0.5 e 
V. Our newly proposed intermediate P^Bj^ differs from the CT state P^BJ^ by a 
compensating charge transfer excitation from Bl to P so that two triplets are created 
which combine to a singlet. Coulomb and exchange energies overcompensate in this 
case the orbital excitation energy to bring the resulting trip-trip singlet P^B]^ far below 
the CT state P+S^and even below P*. The couplings between P^Bj^ to P* and to 
P^Hj^ involve two electron transfer overlaps, which differ from the Dexter exchange 
couplings relevant for excitation energy transfer only in the partitioning of different 
orbitals. The overall tunneling process from P* to P^H£ can be looked upon as an 
electron assisted tunnehng in the sense that one electron from the HOMO of Bl is 
moved to the LUMO of P and back once the transferred electron has passed B^. This 
way the effective tunneling barrier is lowered (Fig.2). 

In order to check the relevance of this mechanism almost all experiments related to the 
charge separation, to recombination and to their magnetic and electric field dependence 
have to be reinterpreted. Particularly informative are changes in the kinetics due to 
mutations or chemical modifications of the pigments. In this paper we shall discuss 
briefiy those cases, which have led to unexpected results within the standard model 
with P^Bj^ as intermediate. Of course the energetics of the trip-trip singlet P^Bj^ is 
sensitive to changes in a different way than that of P^B^. Crudely speaking it is only 
weakly dependent on changes of the electrostatic potentials but more sensitive to local 
distortions of the pigments structure. We will show that these characteristics of our 
model help to interpret the experimental findings including double mutants under a 
new perspective without multiparameter adjustments (Bixon et al 1996;). 

2 Energetics of charge transfer states 

In an attempt to interrelate structure and function for the reaction center it is instruc- 
tive to analyze the energetics of the low lying charge transfer states with regard to the 
following five contributions: 

a) the local structure of the isolated pigments, which is used for the calculation of the 
ionization potentials and the electron affinities, 

b) the Coulomb interactions between the prosthetic groups, which are induced by a 
charge transfer transition, 

c) the polarization effects of these pigments which are due to charge reorganization 
under the influence of the induced Coulomb forces, 

d) the electrostatic polarization effects of those residues and water molecules which are 
in close neighborhood to the prosthetic groups, 

e) the long range electrostatic effects resulting from polar groups of the protein. 

In Fig. 1 those molecules are shown which are treated quantum mechanically. To get 
the local structural effects (a) the ionization potentials (IP) and the electron affinities 
(EA) were evaluated for the six isolated pigments Pm , Pl , Bm , Bl , Hm and Hl- We 



used a semiempirical MO program of the INDO type with the parametrization similar to 
Zerner's ZINDO method (Thompson and Zerner 1990) with configuration interaction 
(Scherer and Fischer 1990) including up to half a million states. The structure of 
Rps. viridis was taken from the protein data bank. The positions of the hydrogen 
atoms were optimized with the help of an MNDO program. The charge transfer (CT) 
induced Coulomb interactions (b) between a donor D and an acceptor A and the polar 
surrounding molecules M have three contributions 

Ec{D+, A~,M)- Ec{D, A, M) = 
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= Ec{D+, A") + Ec{D+, M) + Ec{A', M) 

While the first term stands for the common point monopole interactions of the CT 
induced charges between the donor D and acceptor A only, the other two account for 
the induced interactions of the donor and the acceptor with the other molecules M 
respectively. The energies for the set of CT states Pi^P^, P^BJ^, Pm^l : Pl^Z-, 
PmHl 1 ^-iid those with L and M exchanged, are shown in Fig. 3. The monomer 
contribution, defined as IP-EA, is smallest for the CT state P^H£, indicating that the 
positive charge is better localized on P^ than on Pm- The calculated electron affinity 
of Hl exceeds that of Bj^ by 0.62 eV. This is much more than found in solution (0.3 eV, 
Fajer et al 1975) . To test our calculation we simulated the molecules in solution by 
relaxing their structure and attaching water to Mg instead of the histidine. This way 
we could approximately reproduce the experimental value of the difference in electron 
affinities and the measured absorption spectrum of Bj^ in solution (Scherer 1990). So 
our calculations predict that the electron affinity oi Bl is much smaller in the reaction 
center than in solution. 

The influence of the Coulomb energy Ec (Eq. 1) introduces a sizable contribution to 
the asymmetry of the two symmetry related CT states PlH£ and PmH^ with the 
first falling now below the latter by 0.44 eV. This is mainly due to the difference in 
the anion interaction of H£ or H^^ with the respective ground state dipoles of P^ and 
Pm {Ec{A~,M) from (1)). The polarization Ep was evaluated with a super molecule 
approach, first for the hexamer, consisting just of the pigments Pl, Pm, Bl, Bm , 
Hl and Hm with the histidines attached to the Mg atoms. We then included the 
short range polarization effects Eg resulting from the neighboring residues and water 
molecules. The largest short range effect evolved for the anion H£ via GLU L104. The 
protein electrostatics Ei has been evaluated by means of the DelPhi-program solving 
the Poisson Boltzmann equation (Scherer et al 1995). It lowers the CT state PlH£ so 
far that it falls below the P* state. This way P^H£ becomes energetically accessible 
for the photo-induced charge separation without activation. In fact it is the only CT 



state which fulfills this condition. The state P^Hj^ is higher by 0.4 e V and the state 
P^BJ^ is above P* by 0.7 eV (vertical energy difference). 

The internal CT states P^Pm ^-nd Pm^l always below P^Bj^ . This is easy to 
understand, since they experience a larger Coulomb attraction. For these states the 
Coulomb effect is shown only together with the hexamer polarization. The surround- 
ings do not narrow the gap between P* and P^B^ substantially. TYR M 208 reduces 
it only by 0.11 e V (Scherer et al 1995; Alden et al 1996) and a similar reduction 
was found for the water molecules (Scherer et al 1995). Interestingly the energy of 
P^B^ remained almost unaffected by the long range electrostatic interactions of the 
protein. This conclusion is in hne with the results of Marchi et al (1993). We could not 
reproduce certain results based on dynamics simulations by Warshel et aP , see also 
(Parson ct al 1990; Alden et al 1995). Apparently they allowed for drastic changes in 
the structure. Following our theoretical prediction for the energy of P^B^, we like to 
rule out this state as a real intermediate. Moreover, the electron transfer coupling be- 
tween P* and P'^B2 is so small that a superexchange mechanism cannot be operative 
either. On the other hand, we will see that the trip-trip singlet P^Bj^ can fulfill the 
basic requirements for the energetics and the couplings. 

3 The trip-trip singlet P^B^ 

Double excited singlet states resulting from two locally excited triplets play a major role 
in biological photosystems. They are optically forbidden but can be populated via rapid 
internal conversion processes. In polyenes, they are assigned as Ag states. For larger 
chains their energy falls below that of the lowest optically active Biu state (Hudson and 
Kohlcr 1972; Hudson and Kohler 1973). In biological systems the carotenoids make 
use of the rapid conversion process. A similar Ag state plays an important role in the 
primary isomerization of the retinal chromophore in bactcriorhodopsin (Schultcn ct al 
1995; Takeuchi and Tahara 1997). Another application is found in the isomerization 
of provitamin D (Sobolewski 1994; Fu et al 1997). 

We were recently able to assign a trip-trip singlet within the reaction center, there 
denoted as double triplet P^Pm (Fischer and Scherer 1997). It has been observed in 
the transient excitation spectrum of P* (Wynne et al 1996). It gets intensity mostly 
from the internal CT state P^Pm ■ Its energy location is 0.35 e V above P*. 

Instead of presenting already here detailed quantum calculations on the corresponding 
trip-trip singlet P^Bj^, we like to relate its energy to the experimentally detected trip- 
trip singlet P^Pm in ^ semiempirical way. Within the dimer the triplet energy is lower 
for Pl compared to Pm- This follows from the analysis of ADMR studies (Hoff and 
Vrieze 1996). The triplets of Pm and of i?L are according to our calculations similar 
in energy, so that the sums of the locahzed triplet energies for P^Pm and for P^Bj^ 
respectively, should be about the same. Their induced Coulomb energies, which result 
from changes in the charge distribution due to the trip-trip-singlet excitations, are very 
different. For PmPl found a strong repulsion of 0.57 eV due to the charge shifts 
into the rings I from the rings HI and due to the head to head arrangement of the 



dimer. For Pj^Bj^ we have a head-tail arrangement which resuhs in a weak Coulomb 
attraction of -0.03 e V. The state Pj^B"[ induces an even larger Coulomb attraction 
of - 0.48 eV. This energy gain is overcompensated by the lower local triplet energy of 
Pl for the P^Bj^ state. Since for P^Bj^ the triplet component of P should be largely 
localized on Pl as it is for alone (Hoff and Vrieze 1996) we predict the final energy 
of P^BJ^, which incorporates CI interaction, close to P*. 

In Fig. 4 the calculated transient spectrum of the trip-trip singlet P^Bj^ is shown to- 
gether with the ground state spectrum and experimental results (Dressier et al 1990). 
They show the change in absorbance for the intermediate between P* and P^H£ rel- 
ative to the ground state absorption. The experimental spectrum applies to the two 
step transfer scheme, whereby the second step is the faster. The calculation is based 
on the tetramer Bm Pl Pm Bl with the histidines attached to the Mg-atoms. Exci- 
tations on Bm are not incorporated. They should largely cancel out in the difference 
spectrum. In the ground state spectrum the strong dimer excitation P* is predicted 
at 990 nm in quite good agreement with experiments. The next band polarized almost 
perpendicular (83°) to P* (0°) is the so called upper dimer band. It is somewhat too 
high in energy in comparison with experiments and so is the B^ excitation, which has 
the polarization angle of 33°. The proper shifts to lower energies are indicated by 
horizontal arrows in Fig. 4. 

The transient P^Bj^ absorption shows a transition at 1060 nm, which correlates nicely 
with the observed transition at 1050 nm (Dressier et al 1990; Zinth et al 1996). Also 
its polarization of 30° is in hne with the observed 28°. In our calculation it is composed 
of two accidentally degenerate transitions on the dimer. Both are difficult to assign 
orbital wise, since many CI components contribute. The situation is similar to the 
trip doublet of P^ (Fischer and Scherer 1997) where the charge is replaced by the 
triplet excitation on Pl. The next transition at 915 nm is localized on Bl- It should 
not be as broad as those on P, which have more charge transfer character. Therefore 
it should be worth while to search for such a state in this energy regime to test our 
model. The dominating strong transition at 750 nm is largely the excitation in 
the presence of the triplet on Pl- Together with the transition of Bl the main change 
in the sign of the observed spectrum around 800 nm can be explained once we allow 
for the shift of B^ and the upper dimer band, which are both experimentally justified. 
In addition a shift of the Bm * excitation (not in the calculation) to higher energies 
is expected which might contribute to the amplitudes of the two components around 
820 nm. Finally there are the two almost perpendicular polarized transitions at 680 
nm followed by one localized on Bl at 665 nm with a small polarization angle of 11°. 
This is also consistent with the observation in this frequency regime. The transitions 
at 1050 nm and at 650 nm have been taken as evidence for the occurrence of (Zinth 
et al 1996). We argue this coincidence might be in some parts accidental, since the 
spectrum of Bj^ in solution is not necessarily representative for the spectrum of the 
P^Bj^ state in Rps. viridis. To prove this point, we evaluated the transient absorption 
spectrum of P'^B^ (Fig 5a) for the isolated pigments P and Bl in the presence of the 
counter charge. There is no transition localized on Bl around 1050 nm. Apparently 
the spectrum is shifted relative to the solution spectrum (Fajer et al 1973) to higher 



energies. For a relaxed structure with a water molecule attached we do find the low 
energy transition consistent with earlier calculations (Scherer 1990) in close agreement 
with the solution experiments (Fajer et al 1973) (Fig. 5b). We think a shift by 0.6 eV 
is significant and it is outside the uncertainties of the calculation and the structure. 

For the coupling between P* and P^Bj^ we have only preliminary results. P^Bj^ can 
be approximated by the localized state Bj^ . It couples to the excitonic component 
of P* which is localized on P^ and to the internal charge transfer states of the dimer. 
The dominant matrix element of the electron-electron interaction reads 



We evaluated the relevant two center integrals with atomic Hartree-Fock wave functions 
and obtained for V as one major contribution 5 cm-1, not so different from the one 
particle coupling for the P* — > P'^B^ based on the INDO approximation (10 cm-1, 
Fischer and Scherer 1987; Scherer and Fischer 1989a; Scherer and Fischer 1989b). In (2) 
there are almost no interferences between the dominating atomic contributions. This 
makes the coupling competitive to the one particle coupling of P* to P^Bj^, which 
contains strong interferences. The P^B-[ P~^H£ coupling involves three center 
integrals which we cannot handle at this time sufficiently well. It must be efficient in 
order to assure the fast kinetics for the second step. The second rate process becomes 
in our model the first charge separating step describing the simultaneous shift of two 
electrons from P to Bl and from Bl to Hi respectively. This implies that modifications 
on P can affect this second rate process not only via changes in the oxidation potential 
of P, but also via a change in the coupling. A similar coupling might be responsible 
for the very rapid charge separation following a B^ excitation, which can bypass P* as 
predicted by us (Fischer and Scherer 1987) and recently verified experimentally (Van 
Brederode et al 1997) for a mutant. 

4 Charge separation for mutants and modified re- 
action centers 

As mentioned above the time dependent transient spectra for native reaction centers 
show two clearly separated kinetics for the formation of the charge separated state 
P^H^ . In the spirit of the two step model (Zinth et al 1996) the rate determining 
step P* P^Bj^ is now to be replaced by the P* — * P^Bj^ transition followed by 
the faster step for the P^H£ formation. To test the model it is informative to analyze 
the energetics of modified reaction centers and to test implications for the kinetics. 
Modifications or mutations can change the energetics of the charge transfer states 
P^BJ^ and P^H£ in three ways - via changes of the oxidation potential of P+ or via 
changes of the redox potentials of P^Bj^ or those of P^H£. The corresponding changes 
of the state P^ Bl should be much smaller in all these cases, since energy differences 





between the LUMO's and HOMO's react less sensitively to such modifications than 

the energies themselves. 

Within the first group two double mutants are of particular interest. For L 151 H 
+ M 160 H (Woodbury et al 1994) two additional hydrogen bonds are introduced 
for the dimer, which raises the oxidation potential by 0.14 eV. The double mutant M 
202 H— >L + L 131 L — > H (Laporte et al 1996) introduces a bacteriopheophytin in 
place of Pm and adds a hydrogen bond on the keto group of Pl- It raises the oxidation 
potential by 0.26 eV. For these systems the charge separation process is still operative 
and an intermediate seems to evolve on the 4 ps time scale better detected for the 
first double mutant, which has signatures similar to the transient of P* (Woodbury et 
al 1994). The final charge separation is weakly activated. The common model needs 
to invoke for these systems a super-exchange matrix element for the coupling, which 
should contain a reduction factor of three orders of magnitude in this rate process 
compared to the two step process. There would be no explanation for an intermediate. 
This is not compatible with the observation. Within our model we could allow for a 
small up shift of the energy of the P^Bj^ state and keep it still as a real intermediate. 
Its transient spectrum can become consistent with the observed spectral changes once 
the rapid bypass process (Fischer and Scherer 1987; Scherer and Fischer 1989b; Van 
Brederode 1997) is incorporated for a excitation. 

Reduction of the oxidation potential by 0.08 eV has been achieved for the L168 H-^F 
mutant for Rps viridis (Arlt et al 1996). This change led to a threefold increase of the 
rate determining step, which enhances the population of the intermediate. To us it 
seems difficult to envision such an increase in the P* P^BJ^ transition since that 
coupling is localized in the region of closest approach between P and Bl. Our new 
coupling mechanism becomes sensitive however to vibronically induced couplings of 
the acetyl group of P^, since this modifies the internal CT character of P*. 

The variations in the kinetics caused by strong modifications of the redox potential of 
P'^B'l are also difficult to envision within the common model. The replacement of 5^ 
by a bacteriopheophytin (Zinth et al 1996) should lower the redox potential by at least 
0.3 eV (Fajer et al 1975). We predict an even larger effect due to structural relaxation. 
This value is larger than the full energy difference between P* and P^H£, so that 
P^BJ^ (BPheo) should form a trap. The same should apply for the Ni replacement 
(Haberle et al 1996) of Mg on Bl. A change of the redox potential of 0.29 eV is 
expected on the basis of measurements in solution. Within our model in both cases 
only minor changes of the energy of the P^Bj^ state are predicted which is consistent 
with the observation of an almost identical decay kinetics for P* (Haberle et al 1996). 
For the recombination the state P^B^ (BPheo) might play in our model the same 
role of an equihbrated state as postulated in the literature (Zinth et al 1996), since we 
would predict its energy location close to that of P^H^. 

Heterodimers are mutations which can affect the energetics of the P^Bj^ state in a well 
defined way (McDowell et al 1991). Since the triplet energy of bacteriopheophytin is 
higher than that of bacteriochlorophyll and since the triplet of P^ is localized on Pl ,we 
predict a stronger increase of the energy P^Bj^ for the L heterodimer as compared to 
the M heterodimer. This is also consistent with the observations (McDowell et al 1991). 



The change of the redox potential of P^HJ, by the replacement of Hl by a pheophytin 
for Rb. sphaeroides (Schmidt et al 1995) does not require a new interpretation as long 
as the P^BJ^ state is replaced by our P^B]^ state. In analogy we would conclude that 
the trip-trip-singlet should be about 450 cm~^ below P*. 

For many other mutations such as the 3 vinyl - 13 OH bacteriochlorophyll substitution 
of Bl (Finkele et al 1992; Nagarajan et al 1990; Finkele et al 1990; Shochat et al 1994) 
the variation of the energetics of P^B^ and P^Bj^ should change in the same direction 
even though different couplings are responsible. 

5 Summary 

Our newly introduced intermediate trip-trip singlet P^B^ may bring us to a better 
understanding of the very special features of the reaction center. Within this model, 
we argue that the real trick for the photo-induced charge separation accomplished in 
the evolution process by the construction of the reaction center is the rapid delocal- 
ization of the excitation energy of P* over the dimer and the monomer in form of the 
trip-trip singlet P^Bj^. The formation of the charge separated state P^H~[ avoids the 
appearance of a radical pair between neighboring molecules. This way the recombina- 
tion is suppressed and the Coulomb attraction between the newly born radical pair is 
already strongly shielded by the intermediate Bl- Moreover, this mechanism avoids 
large nuclear reorganization for the initial step, since the state P^Bj^ undergoes no 
strong dipole change with respect to P*. For P^H~[ the Coulomb interaction is suf- 
ficiently shielded. We have shown that the detailed engineering makes largely use of 
the electrostatics, most important of that between the prosthetic groups. The unidi- 
rectionality is partly caused by the surrounding but mostly by the asymmetry within 
the dimer which localizes the triplet and to some extent also the positive charge on P^. 

An additional way to test the model might come from the magnetic field effects. The 
couplings for the triplet recombination are altered in this model and different predic- 
tions result for the Ni mutant (Haberle et al 1996) within the two models. In particular 
the STiperexchange coupling for recombination of the radical pair P^H£ to the triplet 
P^ should be differently affected. 
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Footnotes 

1 A.Warshel kindly provided us with a protein conformation from his dynamics simu- 
lation which put the energy of P^Bj^ below that of P*. In this structure some residues 



were displaced by several A and the Mg-HIS bond at Pl was interrupted by an H- 
atom. 
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Figure Captions 

Figure 1: The quantum chemically treated part of the reaction center Rps. viridis is 
shown. It consists of the four bacteriochlorophylls Bm , Pl, Pm and B^, the two 
bacteriopheophytines Hm and Hl, the quinon Qa and several protein residues 
which are in close contact. 

Figure 2: The electron assisted electron transfer mechanism is visualized. The solid 
arrows refer to the initial presumably slower (3.5 ps) process P* P^Bj^, which 
invokes a two electron exchange between P and Bl. The dashed arrows give the 
second faster (0.65 ps) charge separating step as a simultaneous transfer of an 
electron from B^ to and one from Bl to P. The orbital energies refer to the 
neutral ground state of the hexamer from Fig. 1. 

Figure 3: The calculated energies of the charge transfer states are shown. Several 
contributions are presented separately. lonisation potentials of the donor IP 
and electron affinities of the acceptor EA are estimated from the calculated MO 
energies of the isolated chromophores. The Coulomb energy Ec results from the 
corresponding electron densities of the MO's. The polarization contribution Ep is 



evaluated from a calculation including the six chromophores as a supermolecule. 
Short range interactions Eg refer to the residues shown in Fig.l and are treated 
explicitly whereas long range electrostatic effects Ei are treated in a continuum 
approximation. 

Figure 4: The calculated transient difference spectrum of P^Bj^ (bars) is compared 
with experimental values (circles) from (Dressier et al 1990). The numbers show 
the calculated polarization angles relative to the dimer band P*. The horizontal 
arrows, indicate the energy shifts needed for the experimental assignment. 

Figure 5a: Calculated transient spectra for the CT state P^B^ based on the structure 
of Rps. viridis. The angles shown for the transitions of B^ and P+ are relative 
to the dimer transition P* . The spectra of P"'"(thin bars and diamonds) and 
Bj^ (thick bars and triangles) are calculated in the presence of the corresponding 
counter charge. They are superimposed. 

Figure 5b: The spectrum of P^ (thin bars and diamonds) is superimposed with the 
spectrum calculated for a relaxed BChl~ H2O anion complex (thick bars and 
triangles) . 
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